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Abstract: 

Top quark pair production by (polarized) yy collisions offers an interesting testing ground of the Standard 
Model and its extensions. In this Eetter we present results for differential cross sections of top quark pair 
production and decay including QCD radiative corrections. We take into account the full dependence 
on the top quark spins. We give analytic and numerical results for single and double differential angular 
distributions of tt decay products which are due to top quark polarizations and spin correlations in the 
intermediate state. 
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I. Introduction 


At a future linear lepton collider, backscattered laser light may provide very high-energy photons [1], 
which would allow for a very interesting physics program [2,3]. In particular, top quark pair production 
is possible with large rates in (un)polarized photon photon fusion. The measurement of the process 
TY —> ttX is an important test of the Standard Model (SM). The hrst order QCD corrections to this 
process have already been calculated and found to be large [4-7]. The electroweak virtual plus soft- 
photonic 0(a) corrections are also known [8]. This process will also provide information on possible 
anomalous ytt couplings [9,10] without contributions from Ztt couplings present in e+e collisions. 
Once the Higgs boson is discovered, it will be of primary importance to determine whether its properties 
are as predicted within the SM. In this respect, the process yy —> ttX may play an important role. For 
example, heavy quark production in polarized yy collisions will help to determine the parity of the Higgs 
boson produced as a resonance and decaying into top quark pairs [6,7]. In particular, if a Higgs boson is 
no CP eigenstate, spin correlations of the top quark pairs will help to probe the scalar and pseudoscalar 
couplings of the Higgs boson to the top quark [11,12]. 

For this kind of studies, predictions for top quark pair production and decay at a photon collider must be 
as precise as possible within the SM. In particular, the spin state of the intermediate tt pair must be taken 
into account. (The role of the top quark polarization in probing the tt threshold dynamics in yy collisions 
was discussed in [13]). 

The purpose of this paper is therefore to study the processes 

yy^tfA^f + / + A, (.+ jets + A, all jets, (I.l) 

where (. stands for a charged lepton, with polarized photons from backscattered laser beams. We in¬ 
clude QCD radiative corrections and take into account polarization and spin correlation effects of the 
intermediate tt pairs. 

Leading order results and QCD corrections for the cross section and for top quark spin observables in the 
process yy^ ttX are summarized in sections [n] and [nil Numerical results to order for the effective 
lepton collider cross section and for several decay distributions are given in section |TV| 


II. Kinematics and leading order results 

The production of top quark pairs by photon scattering at leading order is described by the reaction 

y(pi,?ii)-f y(p2,?^2) -^t{kx,st)+t{k2,st). (II.l) 

Here, p\,p 2 ,k\ and k 2 denote the momenta of the particles, Xi and ?i 2 are the helicities of the photons, 
and the vectors St and sj describe the spins of top quark and antiquark. These fulfil the relations 

= Sj = — \ and ki- St= kj ■ sj = 0. (IL2) 

In the (anti)top rest frame the spin of the (anti)top is described by a unit vector s? (Sf). We choose the 
specific rest frames that are obtained by a rotation-free Lorentz boost from the zero momentum frame 
of the tt quarks (tf-ZMF). Both the tf-ZMF and the t and i rest frames will be used to construct spin 
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observables from the final state momenta of the tt deeay products. We use the tf-ZMF rather than the 
c.m. frame of the colliding high-energy photons, since the latter system is probably more difficult to 
reconstruct experimentally. For the 2^2 process of Eq. (III. Ill , the two frames coincide. 

The differential cross section for the process of Eq. dll.ll) can be written as follows: 

d<5(klM,St,St) = ;^|Mopjr 2 , (II.3) 

ZiS’yy 

where the two-particle phase space measure is denoted hy dT 2 , 8^^= and = 3 is the number 

of colours. A simple calculation gives: 

_ [ 32 ^ 2^2 1^0 + -Bo [p\ ■ (■^f + ■5'?)] +5oki^A.2 [P2 ■ (■s'r + ■5'?)] 

-F CQ{sfSt)+DQ{pi-St){p2-St)+DQ\z^-z{p2-St){p\-St)\ (II.4) 


with 


Ao = l+2f{l-z^)-^''[l + {l-zY]+lil2[l-2f{l-z^)-^\^{2-z^)], 
Bo = — [^i(l-2p2 + pV)+X2(l-13V)], 

^yy 

Co = l-2p2 + p4^1 + (l-^2)2]^;,^;,2[l-2p2 + pV(2-z2)], 

4(l-h|3z)(l-z2)(l-XiX2)p2 
L>o =-• 

Sjy 

Here, Qt = 2j3, mis the top quark mass. 


(11.5) 

( 11 . 6 ) 

(11.7) 

(11.8) 



and z is the cosine of the scattering angle in the tt-ZMF, i.e. 
one of the photons (of the top quark) in that frame. 


(11.9) 

Z^Py-k, where Py (k) is the direction of 


III. NLO results for yy ^ ttX 

In this section we present results for the inclusive reaction 

YY^ttX (III.l) 

to order a^a^. Apart from the cross section we study observables that depend on the spins of the top 
quark and antiquark. For polarized photons, observables of the form 


0^ = 2Sta 


( 111 . 2 ) 
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can have non-zero expectation values. Here, a is an arbitrary reference direction and Sf is the top quark 
spin operator. The expectation value of is related to a single spin asymmetry: 


(a') 


g(T)-g(i) 
o(T) + aa)’ 


(III.3) 


where the arrows on the right-hand side refer to the spin state of the top quark with respect to the 
quantization axis a. We will consider here two choices for a, 


a = k (helicity basis), 

a = p (beam basis), (111.4) 

where k denotes the direction of the top quark in the tf-ZMF and p is the direction of the lepton beam 
coming from the left in that frame, which coincides to good approximation with the direction of one 
of the high-energy photons. Top quark polarization perpendicular to the plane spanned by p and k is 
induced by absorptive parts in the one-loop amplitude. This effect is, however, quite small (~ a few 
percent) [14]. 

Analogous observables may of course be defined for the top antiquark. 

Apart from the above single spin observables, we also consider observables of the form 

6 >^ = 4(Sfa)(Sfb). (III.5) 

Here, a and b are arbitrary reference directions and Sj is the top antiquark spin operator. The expectation 
value of is related to a double spin asymmetry: 


^ o(TT)+o(U)-a(U)-aat) 

^ ^ a(TT)+a(ii) + a(U)+a(it)- 


(III.6) 


For the reference directions we will consider here 


a = —b = k, 

a = b = p. (III.7) 


Finally, we also present results for the observable 

= (III.8) 

The above double spin asymmetries {0‘^) and (O^) have also proved useful for an analysis of spin 
correlations of top quark pairs in hadronic collisions [15]. 

The NLO cross section for the reaction yy ttX may be written in terms of two scaling functions: 


o(f,m, Xi,?i2) 



(p,Xi,X2) +47ta,c(^)(p,Xi,X2) 


(III.9) 


Likewise, the unnormalized expectation values of the above spin observables are of the form 





rrP- 


ji®^(p,Xi,X2) +47ta,ji^^(p,Xi,X2)l , 


(III. 10) 
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where a = 1 eorresponds to the observable defined in Eq. (ITTT.8I) . a = 2(3) to the observable O'^ 
defined in Eq. (ITTT-51) in the helieity (beam) basis, and a = 4(5) eorresponds to the single spin observable 
O'’ defined in Eq. (ITTT.2II in the helieity (beam) basis. The variable p is defined as 


AnP' 
P =- 

iyy 


(III.ll) 


The lowest order sealing funetions and ean be eomputed analytieally. We use the following 
auxiliary funetions, whieh vanish in the limit (3 = ^^1 — p —0: 

fi = ^ [ln(ji;) + 2p], 


^2 pa 

= W[ 


ln(x) + 2(3 + -|3^ 




where x = (l — (3)/(l+(3). We then obtain: 

c(®)(p,?ll,?l2) = A^Ttpp jl+p-p^ + ?ll?l2- 


1 + P — ^ — ^1^2 


h 


<)(p,xaz) = 


df(pVl.>-2) = WICPpI 


+ p + p^ + (l + 2p)A.iA.2 + 
l+8p — 7p^ + p2-|-(5 — 3p + p^)A.iA.2 


1 - —-(l + p)?ll?l2 


,2 I P 


1 + p — 2p H—-—h ( 1 H —— ) ^1^2 


4 , 


4 °^(P,Xi,X 2 ) = A^Ttpp 


9 - 20yp - 6p + 14p3/2 + 29p2 + 6p5/2 - 20p3 + 3p^ 

15 


-21 - 20yp + 25p + 14p3/2 - 16p2 + 6p5/2 - 3p3 
+ -- AiA 2 


+ 


15 

_ 1 _ 4p - p2 - 2yp + ^ - 2p3/2 + p5/2 


+ I 1 + 2 y^ + 3p + — I A 1A2 


H-Vpf-h}. 

df^ip.XiM = A'icpp{yr^-|in(x)}(>.,+>.2). 

df\p,h,X2) = -iVltPp|f^(l-6vi; + pV’- 

i-yp 


+ 


2 + 2a/p + 3p — p^/^ — p"^) £2 (Ai — A 2 ). 


(III. 12) 


(III. 13) 
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Figure 1: Left: Scaling functions (dotted), 1) (dash-dotted), (full), 

and c^^\p,l, — \) (dashed). Right: Scaling functions df\p,\,\) (dotted), j|^^(p,l,—1) (dash-dotted), 
j|^^(p, 1 , 1 ) (full), and ^i?|^^(p, 1 ,- 1 ) (dashed). 




The functions and d^^ are obtained by a numerical integration. The scaling functions for the spin- 
averaged cross section and all spin observables are plotted in Figs. 1-3 for different choices of the photon 
helicities as a function of 

ri = --l. (III.14) 

P 


The result for unpolarized photons can be inferred from 


1 r 


c"»'('>(p,0,0) = r c(»M‘)(p,l,l) + c(»M»(p,l,-l) 


d|“M‘>(p,0,0) 

“ 4,5 


1 L(o),(i) 


(p,0,0) = 


2 

0 . 


(P,1,-1) 


(III. 15) 
(III. 16) 
(III. 17) 


As a check we compared our result for the functions c*^®)’*^^)(p,0,0) with the results given in Fig. 2 of 
Ref. [4] and found perfect agreement. We further compared the functions 1, ±1) to the results 

given in Table 1 of Ref. [ 6 ]. After a trivial rescaling to account for the different conventions used in the 
definition of the scaling functions, we also found agreement. 
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r|=1/p-1 ri=1/p-1 


Figure 2: Left: Scaling functions d^\p, 1,1) (dotted), d^\p, 1,-1) (dash-dotted), 1,1) (full), 

and 1,-1) (dashed). Right: Scaling functions d^\p, 1,1) (dotted), df\p, 1, —1) (dash-dotted), 

J 3 ^^(p, 1,1) (full), and ^i? 3 ^^(p, 1,-1) (dashed). 
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Figure 3: Left: Scaling functions 1,1) (dotted), 1,-1) (dash-dotted), ^i? 4 ^^(p, 1,1) (full), 

and d^\p, 1,-1) (dashed). Right: Scaling functions df\p, 1,1) (dotted), d^^\p, 1, —1) (dash-dotted), 
J 5 ^^(p, 1,1) (full), and ^i? 5 ^^(p, 1,-1) (dashed). 
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IV. Effective cross sections and spin observables 

IV. 1. The effective cross section for 77 ^ ttX 

The total tt cross section at a photon collider may be written at NLO QCD as (cf., e.g., [5]) 




OtT = 


f }^max 


Jo 


dyi J^'^^^dy 2 f^iyi,Pe,PL)fyiy 2 ,Pe,PL) +47ta,c(i)}. 


(IV.l) 


The function is the normalized energy spectrum of the photons resulting from Compton 

backscattering of laser light off the high energy electron beam. It is explicitly given by: 


f^{y,Pe,PL)=X-' 


1 


1-y 


-y+i2r-lf-PePLxr{2r-l)i2-y) 


(IV.2) 


Here, PeiPt) is the polarization of the electron (laser) beam, and y is the fraction of the electron energy 
in the c.m. frame transferred to the photon. It takes values in the range 

0 < y < = ymax, (IV.3) 

x+\ 


with 


ml 


(IV.4) 


where Et{Ee) is the energy of the laser (electron) beam and mg is the electron mass. In order to avoid the 
creation of an pair from the backscattered laser beam and the low energy laser beam, the maximal 
value for x is 


'^^max — 2(1 + V^). 

For a beam energy Eg = 250 GeV, this leads to an optimal laser energy 

El ~ 1-26 eV, 

which will be used in the following numerical results. Finally, 

_ 3^ 

' 41-3^)' 

The normalization factor in Eq. (IIV.2I) is determined by 


(IV.5) 

(IV.6) 


(IV.7) 


rymax 

f^{y,Pg,PL)dy=l. (IV.8) 

The scaling functions have to be evaluated at p = 4m^/ {yiy 2 See) and for polarizations 


l.=Pyiyi,pl‘\pl;^), i=l,2. 


(IV.9) 
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Table 1: Results for the effective cross section at JJfe — 500 GeV. 


{Pel,Pe2',PLl,PL2) 

[fb] 


K = 

(0,0;0,0) 

49.81 

76.44 

1.53 

(0.85,0.85;-!,-1) 

175.86 

260.77 

1.48 

(0.85,0.85; + l, + l) 

15.96 

26.89 

1.68 

(0.85,-0.85;-l, + l) 

48.99 

71.93 

1.47 


The funetion Py{y,Pe,PL) describes the degree of polarization of photons scattered with energy fraction 
y, which is given by 


Pjiy.Pe.Pt) 


1 

^ffiy.Pe.Pt) 


i^xrPe [1 + (1 -y)(2r- 1)2] - (2r- 1)Pl 


1 


+ l-y 


(IV. 10) 


Numerical results for are given in Table [T] for = 500 GeV and different polarizations of the 
laser and electron beam. We use the values = 178 GeV, a = 1/128 and as{iJ = mt) =0.1. The QCD 
corrections to a»- are quite large. This is because for = 500 GeV most of the top quark pairs are 
produced close to threshold where the Coulombic |3 ^ singularity from soft gluons is important. 


IV.2. Spin observables 

The spin observables for yy —ttX discussed in section [nil translate into observables built from the 
momenta of the tt decay products. 

The single spin asymmetries (ITTT.2II cause a nontrivial one-particle inclusive decay distribution of the 
form 


1 ai -hV) 

o JcosGi 



-fB/cosGi). 


(IV. 11) 


Here, Gi is the angle between the direction of a top quark decay product ai measured in the top quark rest 
frame and one of the reference directions a defined in Eq. (ITTT.4II . The coefficient B,, with i = heli,beam 
for the helicity and beam bases, is determined by the top quark spin asymmetry (ITTT.2II and by the so- 
called spin analysing power of the decay product ai, which will be discussed below. 

The double spin asymmetries (ITTT.5II lead to a two-particle inclusive decay distribution of the following 
form: 


Ja(YY 0102+X) -j-B/cosGi -hB,cosG2-C,cosGicosG2) , (IV. 12) 

a acosGiacosG 2 4 

where Gi is defined as above and G 2 is analogously the angle between one of the top antiquark decay 
products and one of the reference directions b defined in Eq. (1111.71) . The coefficients C, are determined 
by the double spin asymmetries (IIII.5I) and the spin analysing powers of the two decay products oi and 
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a 2 . Finally, a non-zero expectation value of the observable defined in Eq. (ITTT.8II leads to a distribution 
of the form 


1 Ja(yy-^ aia2+X) 
a Jcoscp 


1 

2 


(1 — Dcostp), 


(IV. 13) 


where cp is the angle between the direction of flight of the top decay product ai and the antitop decay 
product a 2 defined in the t and t rest frames, respectively. We recall that these rest frames have to be 
obtained by a rotation-free boost from the tf-ZMF. 

The spin analysing power of the t and t decay products is encoded in the one-particle inclusive angular 
distributions dT/d cos 0 for the decays 


t{st) ai{qi) +Xi , 

t{st)^a2{q2)+X2. (IV.14) 


Here qi and q 2 are the momenta of ai and a 2 , respectively, defined in the rest frame of the (anti)top 
quark, and 0 is the angle between the polarization vector of the (anti)top quark and the direction of flight 
of ai (a 2 ). For a fully polarized ensemble of top quarks (antiquarks) these distributions are of the form 


dcosd 


r(h2) 

2 


(1 COS0) , 


(IV.15) 


where is the partial width of the respective decay channel. The quantity is the (anti)top- 
spin analysing power of ai^ 2 - For the case of the standard {V — A) charged current interactions these 
distributions were computed to order for the semileptonic and non-leptonic channels in Refs. [16] 
and [17], respectively. 

As we work to lowest order in the electroweak couplings, F*^^) = F^^) and K 2 = Ki to all orders in a^, if 
the channel a 2 +X 2 is the charge-conjugate of ai -I-Ai. 

For semileptonic top decays t —> the charged lepton is the most efficient analyser of the spin 

of the top quark. In the case of non-leptonic decays t bqq'{g) a good top spin analyser that can be 
identified easily is the least-energetic light quark jet. 

In Ref. [17] the coefficients were given to NFO accuracy for different choices of the spin analyser. 
To compute the coefficients B;, C,- and D we need the partial widths 


p(sl,h) 


_ „(sfh) 

— Uq 


+ AnUsa 


(sl,h) 

1 ’ 


(IV. 16) 


where the indices si and h stand for semileptonic and hadronic decay modes. Further, we need the 
dimensionful coefficients 

p(sl,h)^(£,y) ^ ^ ^ (IV. 17) 

where £ (j) refers to using the charged lepton (least-energetic light quark jet) as spin analyser. For the 
determination of these coefficients we use the Fermi constant Gf = 1.16639 x 10-5 GeV-2, m = 178 
GeV, niw = 80.42 GeV, Fvk = 2.12 GeV, = 4.75 GeV, and all other quark and lepton masses are put to 
zero. (We do not use the narrow width approximation for the intermediate W boson.) We obtain, putting 
the CKM matrix elements \ Vtb\ = \ Vqq'\ = 1: 

al = 0.52221 GeV, 
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a\ = -0.01968(15) GeV, 
af = -0.01097(5) GeV. 


For the relevant eoeffieients ho.i we obtain: 


bl = 0.26950 GeV, 
b\ = -0.01118(8) GeV, 
b{ = -0.02375(26) GeV. 


(IV. 18) 


(IV. 19) 


The Durham algorithm was used as jet clustering scheme to obtain the four parton contribution to b{. 
Within the leading pole approximation for the intermediate top quarks and antiquarks, the coefficients 
of the single and double differential distributions dlV.l 1I) - (IIV.13I) are obtained in terms of the following 
quantities: 


a. = 


Od = 


N. = 


V" = 



aJQf 1 


m2 Ft. 

X 



a?Q^ 1 


m2 r2 

X 



a^Qf 1 


m2 Ft. 

X 



cPQj 1 


m2 r2 

X 



dyi 
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^Vmax /*ymax 

dy2f^(yi,Pe,PL]f^(y2,Pe,PL] 

r (Vmax 

dy\ / dy 2 fy{yuPe,PL)fy{y 2 ,Pe 2 PL) 


+ 47taj 

r rmax 

Jo 
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T “0 


-tQ ^ 


^Vmax 


C ymax 


dyi / dy 2 ayuPe,PL)fJ{y 2 ,Pe,PL) 


+ 47tai 

P (Vmax 


y 


to 


dy\ 


C (Vmax 


dy2myuPe2PL)my22Pe2PL) 


|. 


We then get to NLO in a^: 


D = 


Vf 


nJ 

Cheli = -, 

<3d 


_ ^3 

Cbeam — ; 

<3d 


K 

Bheli = -, 


Bbeam — 




(IV.20) 


(IV.21) 


(IV.22) 


(IV.23) 


(IV.24) 


The LO and NLO results for these quantities are shown in Table |2l for favorable and realistic choices 
of electron and laser polarizations, using the same parameters as in Table [B In most cases the QCD 
corrections are of the order of a few percent and thus much smaller than the corrections to the total tt 
cross section. This was to be expected, since the bulk of the corrections is due to soft gluons which do 
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not affect the tt spin state. The biggest correction (~ 11%) occurs for the coefficient Bbeam if the least 
energetic light jet is used as spin analyser. 

Photon polarization is an important asset: It is necessary to obtain polarized top quarks and thus non¬ 
zero coefficients Bheii and Bbeam- Further, the choice {Pe]_,Pe 2 \PLi,PL 2 ) = (0.85,0.85; —1,-1), which 
increases the total yield of tt pairs by more than a factor of 3 (see Table 1), in addition leads to larger tt 
spin correlations. In particular, in the helicity basis the correlation is then almost 100% in the dilepton 
channel. 


Table 2: Results for double and single spin asymmetries at JJfe = 500 GeV. 



{Pe\,Pe2',PLl,PL2) 

dilepton 

LO NLO 

lepton-jet 

LO NLO 

jet-jet 

LO NLO 

D 

(0,0;0,0) 

-0.670 

-0.686 

-0.346 

-0.338 

-0.178 

-0.167 


(0.85,0.85;-!,-1) 

-0.806 

-0.801 

-0.416 

-0.394 

-0.215 

-0.194 

Cheli 

(0,0;0,0) 

0.811 

0.826 

0.418 

0.408 

0.216 

0.201 


(0.85,0.85;-!,-1) 

0.985 

0.981 

0.508 

0.483 

0.262 

0.238 

Cbeam 

(0,0;0,0) 

-0.580 

-0.606 

-0.299 

-0.299 

-0.154 

-0.148 


(0.85,0.85;-!,-1) 

-0.808 

-0.804 

-0.417 

-0.396 

-0.215 

-0.195 



lepton-i-X 

jet-tX 




{Pe\,Pe2',PLl,PL2) 

LO 

NLO 

LO 

NLO 



Bheli 

(0.85,0.85;-!,-1) 

0.658 

0.655 

0.340 

0.323 



Bbeam 

(0.85,-0.85;-!, 1) 

-0.684 

-0.637 

-0.353 

-0.314 




So-called non-factorizable corrections do neither contribute at NLO QCD to nor to the angular cor¬ 
relations considered above. A proof of this statement is given in [15]. 


V. Conclusions 

We have computed a variety of spin observables for the process ttX up to order a^a^. Together with 
the differential rates of polarized top and antitop quark decays at order a^, we have obtained the NLO 
QCD contributions to the fully differential cross section with intermediate top quark pair production at a 
photon collider. 

We have applied the above results to tt production and decay at a future linear collider operating at 
= 500 GeV. We have shown that for an appropriate choice of the polarizations of the laser and 
electron beam, the cross section and the double/single spin asymmetries can be quite large. While the 
QCD corrections to the cross section can be very large, most of the double/single spin asymmetries 
are affected at the level of only a few percent. The observables considered here will provide useful 
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tools to analyse in detail the top quark pair production and decay dynamics. In particular, their precise 
measurement will test whether the top quark truly behaves as a quasi-free fermion as predicted in the 
Standard Model. 
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